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Abstract 

The rare decay —> Df f can occur only via annihilation type diagrams in the standard model. 
We calculate this decay in perturbative QCD approach with Sudakov resummation. We found that 
the branching ratio of B'^ —> Df f is of order 10“^ which may be measured in the near future by 
KEK and SLAC B factories. The small branching ratio predicted in standard model makes this 
channel sensitive to new physics contributions. 
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I. INTRODUCTION 


The rare B decays are useful for the test of standard model (SM). They are sensitive to 
new physics contributions, since their branching ratios in SM are small. Some of them have 
already been measured by CLEO and B factories in KEK and SLAG. Most of them are 
still on the way of study by both experimental and theoretical studies. Among them, the 
inclusive or semi-inclusive decays are clean in theory, but with more uncertainty in exper¬ 
imental study. On the other hand, the exclusive decays are difficult for precise theoretical 
prediction but easier for experimental measurement. The study of exclusive rare B decays 
require the knowledge of hadronization, which is non-perturbative. The generalized factor¬ 
ization approach has been applied to the theoretical treatment of non-leptonic B decays for 
some years |]l|. It is a great success in explaining many decay branching ratios P, |^. The 
factorization approach is a rather simple method. Some efforts have been made to improve 
their theoretical application IQ and to understand the reason why the factorization approach 
has gone well P). One of these method is the perturbative QCD approach (PQCD), where 
we can calculate the annihilation diagrams as well as the factorizable and non-factorizable 
diagrams. 

The rare decay B^ —> Df(j) is pure annihilation type decay. The four valence quarks in 
the final states Dg and 0 are different from the ones in B meson, i.e. there is no spectator 
quark in this decay. In the usual factorization approach, this decay picture is described as b 
and u quark in B meson annihilating into vacuum and the Dg and 0 meson produced from 
vacuum then afterwards. To calculate this decay in the factorization approach, one needs 
the Dg ^ (f) form factor at very large time like momentum transfer 0{Mb)- However the 
form factor at such a large momentum transfer is not known in factorization approach. This 
makes the factorization approach calculation of these decays unreliable. 

In this paper, we will try to use the PQCD approach to calculate this decay. The W 
boson exchange causes the four quark operator transition bu —>■ sc, the additional ss quarks 
included in DgCp are produced from a gluon. This gluon attaches to any one of the quarks 
participating in W boson exchange. This is shown in Figure |^. In the rest frame of B meson, 
both s and s quarks included in Dg(j) have 0{Mb) momenta, and the gluon producing them 
also has momentum ~ This is a hard gluon. One can perturbatively treat the 

process where the four quark operator exchanges a hard gluon with ss quark pair. It is just 
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the picture of perturbative QCD approach [^, |^. Furthermore, the hnal state of this decay is 
isospin singlet. It is proportional to the Vub transition. Any Vcb transition can not contribute 
to it. Therefore there will not be any dominant soft hnal state interaction contribution. 
Unlike the B ^ KK decays (which may have large contribution from hnal state interaction 
contribution) 0 , the decay B —> Dgip is a very clean channel for a test of annihilation type 
contribution. 


In the next section, we will show the framework of PQCD briehy. In section |T|, we 
give the analytic formulas for the decay amplitude of B^ decays. In section 

we give the numerical results of branching ratio from the analytic formulas and discuss the 
theoretical errors. Finally, we conclude this study in section 0. 


II. FRAMEWORK 


PQCD approach has been developed and applied in non-leptonic B meson decays |^, 
0. B ft 0.13 for some time. In this approach, the decay amplitude is described by three 
scale dynamics, soft(<I)), hard(Ff), and harder(C) dynamics. It is conceptually written as 
the convolution. 


Amplitude 



kid^k2d^h TT[C{t)^Biki)^Dsik2)^4,ik3)H{ki, k2, k3,t)], 


( 1 ) 


where /cj’s are momenta of light quarks included in each mesons, and Tr denotes the trace 
over Dirac and color indices. C{t) is Wilson coefficient of the four quark operator with 
the QCD radiative corrections. C{t) includes the harder dynamics at larger scale than Mb 
scale and describes the evolution of local 4-Fermi operators from W boson mass, down 
to f ~ 0{Mb) scale, which results in large logarithms \n{Mw/t)- H describes the four 
quark operator and the quark pair from sea connected by a hard gluon whose scale is at the 
order of Mb, and includes the 0{Mb) hard dynamics. Therefore, this hard part H can be 
perturbatively calculated, t is chosen as the largest energy scale in H, in order to lower the 
corrections to hard part H. is the wave function which describes hadronization of 
the quark and anti-quark into the meson M. While H depends on the processes considered, 
is independent on the specific processes. Determining in some other decays, we can 
make quantitative predictions here. 
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We consider the B meson at rest for simplicity. It is convenient to use light-cone coordinate 
iP~ iVt) fo describe the meson’s momenta, where = (p° ±p^)/-\/2 and pr = 

On this coordinate we can take the B, Dg, and 0 mesons momenta as Pi = "^(1,1, Oi^), 
F 2 = ^(l,r 2 , 0 r), and P 3 = ^( 0,1 - r^Or), respectively, where r = MdJMb and we 
neglect the square terms of 0 meson’s mass M|. Putting the light spectator quark momenta 
for F, Dg and 0 mesons as /ci, ^ 2 , and /C 3 , respectively, we can choose ki = ( 0 , xiPf, kir), 
^2 = 0, k 2 T) and ^3 = {0,X3P^,'k 3 T)- Then, integration over , k^ and k^ in eq.(|l]) 

leads to 


Amplitude ~ Jdxidx 2 dx 3 bidbib 2 db 2 bsdbs 

Tr [0(t)<I)B(xi, 6 i)$D,(a; 2 , & 2 )*h^(a; 3 , bii)H{xi, h, , (2) 

where bi is the conjugate space coordinate of ki^. The last term, e~^, contains two kinds 
of logarithms. One of the large logarithms is due to the renormalization of ultra-violet 
divergence In tb, which describes the QCD running between scale t and 1/6. The other is from 
double logarithm due to soft gluon corrections. This double logarithm called Sudakov form 


factor suppresses the soft dynamics effectively [^. Thus it makes perturbative calculation 
of the hard part H applicable at intermediate scale, i.e., Mb scale. We calculate the H for 
P+ ^ D^(j) decay in the hrst order in Ug expansion and give the convoluted amplitudes in 
next section. 

In order to calculate analytic formulas of the decay amplitude, we use the wave functions 
"I’m, 0/3 decomposed in terms of spin structure. As a heavy meson, B meson wave function is 
not well dehned. It is also pointed out by the recent discussion of B meson wave function 


131 that, there is no constraint on B meson wave function, if three-parton wave functions 
are considered. To be consistent with previous calculations i,|, pd]| , we follow the same 


argument that the structure ( 7 ^ 75 ) 0/3 and 750 ,^ components make the dominant contribution 
in B meson wave function. Then, <hM,o /3 is written by 


$ 


M,al3 




{(Z’ m75)q:/ 30M + 75a/30M} ) 


(3) 


where A/ = 3 is color’s degree of freedom. Pm is the corresponding meson’s momentum, and 
0^^ are Lorentz scalar wave functions. As heavy quark effective theory leads to 0^ ~ Mb0b, 
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then B meson’s wave function can be expressed by 


$s,a/3(x, h) = -J= [fi + Mb] 75a/30B(a;, b). 


( 4 ) 


According to Ref.|[14||, a pseudo-scalar meson moving fast is parameterized by Lorentz 
scalar wave functions, 0, (pp, and (pa- as 

rl 


{Ds (^)k(^) 7 M 75 c( 0 )| 0 ) ~ -ifDsPf. / dx Xx), 


{P>s (^)k(^) 75 c( 0 )| 0 ) = -ifos^oDs / dx "(pp{x), 


(5) 

( 6 ) 


{Ds iP)\s{zh5(^puci0)\0) = -fn.moDs ( 1 - 


Ml 


m, 


ODs 


{PpZy - PyZp) i dx P'^^^Cpai.x), (7) 


where moD^ = Mj^Jimc + rris)- We ignore the difference between c quark’s mass and Dg 
meson’s mass in the perturbative calculation. This means, = rriQBsi this approxima¬ 
tion, the contributions of eq.(0) are negligible. With the equation of motion eq.(|^), eq.(^), 
lead to 


(pp{x) = (p{x) + O 


A 


M, 




( 8 ) 


Therefore the Dg meson’s wave function can be expressed by one Lorentz scalar wave func¬ 
tion. 


h) = —= [(75 f2)a/3 + MD,-f5a^] b). 

V 

The wave function for M = B, Dg meson is normalized by its decay constant Jm 

[ dx ^m{x, b = 0) = 


(9) 




( 10 ) 


In contrast to the B and Dg meson, for the (p meson, being light, the component 
remains. In B~^ decay, the (p meson is longitudinally polarized. Then, cp meson’s 

wave function is parameterized by three Lorentz structures 




4=4>*fe). 






(11) 


In the numerical analysis we will use $<;(,, and which were calculated from QCD sum 
rule 1^. They will be shown in section |^. 
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FIG. 1: Diagrams for —> D^(f) decay. The factorizable diagrams (a),(b) contribute to Fa, and 

non-factorizable (c), (d) do to Ma- 

III. PERTURBATIVE CALCULATIONS 


The effective Hamiltonian related to D '^0 decay is given as 


15 


O 


+ C2{p)02{p)] 

1 = (lri^PL!l)(n'‘PLu), O 2 = (67pP[,ti)(c7'‘FLs). 


( 12 ) 

(13) 


where C'i^ 2 (/^) are Wilson coefficients at renormalization scale fi. The projection operator 
is defined as Pl = 1 — 75 . The lowest order diagrams contributing to are drawn 

in Fig.|l| according to this effective Hamiltonian. As stated above, B^ —> D'^(f) decay only 
has annihilation diagrams. 

We get the following analytic formulas by calculating the hard part H at hrst order in 
as- Together with the meson wave functions, the amplitude for the factorizable annihilation 
diagram in Fig.|^(a) and (b) results in, 

pi poo P 

Fa = IQhCfIeMI / dx2dx3 / b2db2hdb^ ^Ds{.X2,b2) x {a;3<F^(a;3, 63) 


r{2x^ - l)r0<Fji(a;3,63) + r(l + 2a;3)r^<F^(x3, 63)}P/(t^)h„(a;2,0:3, 62, ^3) 

- {x 2 ^ 4 >{xz, bs) + 2 r(l + X2)r^^l{x3, b3)}Ef{tl)ha{x3, X 2 , b^, 62 ) 


(14) 


where Cp = 4/3 is the group factor of SU(3)c gauge group, and The function 

Ef, ha are given in the Appendix. Since we only include twist 2 and twist 3 contribu¬ 
tions in our PQCD approach, all the and terms in the calculation are neglected for 
consistence. The explicit form for the wave functions, ^m, is given in the next section. From 
eq. (0), one can see that the factorizable contribution Fa is independent of the B meson wave 
function, but proportional to the B meson decay constant fs- 
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The amplitude for the non-factorizable annihilation diagram in Fig.|^(c) and (d) is given 


as 


Ma = / dxidx2dx3 / bidbib2db2 &i)<hn),(x 2 , 62 ) 

Jo Jo 


X 


{x2<h0(a:3, 62) +r{x2- X3) r^$*^(x3, 62) 


+ r (0:2 + 2:3) r^^^(x3, b2)}Bm(tj^)hi^'> (xi, X2, X3, bi, 62) 

- {a;3<h^(a;3, 62) - r (0:2 - 0:3) r^^!p(x3, 62) 

+ r (2 + X 2 + X3) r 0 <h^(a; 3 , 62 ) (^ 1 , ^ 2 , X3,bi, 62 ) 


(15) 


Unlike the factorizable contribntion Fa, the non-factorizable annihilation diagram involve all 
three meson wave fnnctions. 

Thus, the total decay amplitude A and decay width T for cf) decay are given as 


A — Fa + Ma, (16) 

i\/f3 

r(B+ ^ D:^) = -^\v:,V^A\^ (17) 

where the overall factor is included in the decay width with the kinematics factor. 

The decay amplitude for CP conjugated mode, B~ —>■ D~(j), is the same expression as 
> D^cf), just replacing V*^Vcs with VubV*^. Since there is only one kind of CKM phase 
involved in the decay, there is no CP violation in the standard model for this channel. We 
thus have Br{B'^ Dt^) — Br{B~ —> D~(p). 


IV. NUMERICAL RESULTS 


In this section we show numerical results obtained from the previous formulas. At the 
beginning, we give the branching ratios predicted from the same parameters and wave func¬ 
tions that are adopted in other works. Secondly, we discuss the theoretical errors due to 
uncertainty of some parameters. 

For the B meson’s wave function, there is a sharp peak at the small x region, we use 


b) = Nbx^{ 1 — xY exp 


X"^ 1 r 

- s'"*")' 


(18) 


which is adopted in ref. i, i pd| . This choice of B meson’s wave function is almost a best 
£t from the B —>■ Kn, tttt. Tip and tiuj decays. For the Dg meson’s wave function, we assume 
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the form as the following, leaving a free parameter 






- a;){l + aD,{l - 2a;)} exp 


-^iu;nbr 


(19) 


This is a rather flat distribution function. Since c quark is heavier than s quark, this function 
is peaked at c quark side, i.e. small x region. The wave functions of the (j) meson are derived 
by QCD sum rules |T^ 


+ 0,21 (3 - + 35^) + 0.69 (^1 + On 

ir ° r 

= 


3^ (4.5 - 11.2a: + 11.2 x^) + 1.38 In ■ 


X 


4v2iVc ^ ' - / ■ 1 - X 

where .^ = 1 — 2x. In addition, we use the following input parameters: 


( 20 ) 

( 21 ) 

( 22 ) 


Mb = 5.279 GeV, Mb. = 1.969 GeV.m* = 1.02 GeV, (23) 

/b = 190 MeV, U = 237 MeV, /J = 220 MeV, fo. = 241 MeV, (24) 

ut = 0.4 GeV, an. = 0.3, ud = 0.2GeV. (25) 


With these values and eq. (pT]) we get the normalization factor Nb = 91.745 GeV. Using the 
above hxed parameters, we hnd that the factorizable annihilation diagram contribution is 
dominant over the non-factorizable contribution. The reason is that the Wilson coefficient 
in non-factorizable contribution is C'i(t), which is smaller than the one in factorizable 
contribution Fa, ai = Gi/d -|- G 2 . Although the real part of Ma is negligible, the imaginary 
part of Ma is comparable with the imaginary part of Fa, it is about 30% of the real part of 
Fa. 

The propagators of inner quark and gluon in Figure 1 are usually proportional to 1/xj. 
One may suspect that these amplitudes are enhanced by the endpoint singularity around 
Xj ~ 0. This can be explicitly found in eq.(^, [A9|) , where the Bessel function Yq diverges 
at Xj ~ 0 or 1. However this is not the truth in our calculation. First we introduce the 
transverse momentum of quark, such that the propagators become l/(xjXj -|- fc^). There is 
no divergent at endpoint region. Secondly, the Sudakov form factor Exp [—S'] suppresses the 
region of small k^. Therefore there is no singularity in our calculation. We also include the 
threshold resummation in our calculation of factorizable diagrams, which further suppress 

















the endpoint region contribution |]^. The dominant contribution is not from the endpoint 
region of the wave function. As a prove, in our numerical calculations, for example, an 
expectation value of in the integration for Fa and Ma results in {ag/ii) — 0.1. Therefore, 
the perturbative calculations are self-consistent. 

Now we can calculate the branching ratio according to eqs. Here we use CKM 


matrix elements 18 


iKbl = 0.0036 ±0.0010, |14,| = 0.9891 ±0.016, (26) 

and the life time for meson is tb± = 1.65 x 10“^^ s. The predicted branching ratio is 

Br(H+ ^ T)+0) = 3.0 x 10“^. (27) 


This is still far from the current experimental upper limit []^ 

Bt{B+ D+(j)) < 3.2 X 10“^ (28) 

The branching ratios obtained from the analytic formulas may be sensitive to various 
parameters, such as parameters in eqs.(p5|). Uncertainty of the predictions on PQCD is 
mainly due to the meson wave functions. Therefore it is important to give the limits of 
the branching ratio when we choose the parameters to appropriate extent. The appropriate 
extent of Ub can be obtained from calculation of semi-leptonic decays fl^ and other B tttt, 
B —>• Kn and B —> pvr, oott decays BBQ. 

0.35 GeV <wt< 0.45 GeV. (29) 

The change of value of will not alter the result of Fa, which is independent of B meson 
wave function, but will affect the value of Ma. We did not hnd any strict constraints for 
the Dg meson wave function in the literature. In fact, a future study of H —> Dgir will 
do this job. At present, in Dg meson wave function is a free parameter, and we take 
0 < uDs < 1- Here we check the sensitivity of our predictions on oob and within the 
ranges stated above. The branching ratios normalized by the decay constants and the CKM 
matrix elements can result in 


Br(H+ ^ Dt<P) = iS.Otli) X 10 


-7 


fs fos 


\v:b Va. 


(30) 


,190 MeV-241 MeVy VO-0036 ■ 0.9891^ 
Considering the uncertainty of fs, fos and etc., the branching ratio of B^ —> D~^(p 

decay is at the order of 10“^. This may be measured by the current B factory experiments 
in KEK and SLAC. 
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V. CONCLUSION 


In two-body hadronic B meson decays, the hnal state mesons are moving very fast, since 
each of them carry more than 2 GeV energy. There is not enough time for them to exchange 
soft gluons. The soft hnal state interaction is not important in the two-body B decays. This 
is consistent with the argument based on color-transparency |^. The PQCD with Sudakov 
form factor is a self-consistent approach to describe the two-body B meson decays. Although 
the annihilation diagrams are suppressed comparing to other spectator diagrams, but their 
contributions are not negligible in PQCD approach Q. 

In this paper, we calculate the B^ —> decay in the PQCD approach. Since neither 

of the bottom quark nor the up quark in the initial B meson appeared in the hnal mesons, 
this process occurs purely via annihilation type diagrams. It is a charm quark (not an 
anti-charm quark) in the hnal states, therefore the usual Vcb transition does not contribute 
to this process. The hnal states are isospin singlet. There should be no dominant hnal 
state interactions through other channels contribute. From our PQCD study, the branching 
ratio of B^ decay is still sizable with a branching ratio around 10“^, which may 

be measured in the current running B factories Belle, BABAR or in LHC-B in the future. 
This may be one of the channels to be measured in B decays via annihilation type diagram. 
Whether the PQCD predicted branching ratio is good enough to account for the B^ —> Dp 0 
decay will soon be tested in the current or future experiments. 

The small branching ratio (comparing to the already measured other B decays) predicted 
in the SM, makes this channel sensitive to any new physics contributions. Since the CP asym¬ 
metry predicted for this channel in SM is zero, any non-zero measurement of CP asymmetry 
will be a dehnite signal of new physics. We also notice that the supersymmetric contribution 
will not enhance the decay branching ratio signihcantly, but it may contribute to a non 
zero CP asymmetry in this channel, since the supersymmetry couplings can introduce new 
phases. 
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APPENDIX A: SOME FUNCTIONS 


The dehnitions of some functions used in the text are presented in this appendix. In the 
numerical analysis we use one loop expression for strong coupling constant, 

Jtt 




(Al) 


/?olog(/xVA 2 )’ 

where Po = (33 — 2n/)/3 and Uf is number of active quark flavor at appropriate scale. A is 
QCD scale, which we use as 250 MeV at n/ = 4. We also use leading logarithms expressions 


for Wilson coefficients Ci ^2 presented in ref. [^. Then, we put m* = 170 GeV, mw = 80.2 
GeV, rrih = 4.8 GeV, and rrtc = 1.3 GeV in the Wilson coefficients calculation. 

The function Ef and are defined as 


Ef{t) = [C'i(f)/3 + C'2(f)]a.(f)e-^-W-5At)^ 


(A 2 ) 

(A3) 


The above are defined as 


SbP) = s{xiP^,bi) + 2 j 

Jl/bi h 

Snit) = s{x2P2,b3) + 2 j %7g(h), 

Jl/b2 h 

S^{t) = S{X3P3, bs) + s((l - X 3 )P;^, 63 ) + 2 [ ), 

Jl /63 


(A4) 

(A5) 

(A 6 ) 


where the last terms of the above formulas are logarithms from the renormalization of ultra¬ 
violet divergence. The term s{Q,b), the so-called Sudakov factor, result from summing up 
double logarithms caused by collinear divergence and soft divergence. The expression is 
given as 0 

'l/b , 


s{Q,b) = [ ^ |^( 27 E-l-log 2 ) + Gi.log^| 
ii/feh L13 h J 


67 7r2 10 2 7 e 

9 3 27^^ ^ 3^° 2 


TT 

71 J p' 


(AT) 


-Ye = 0.57722 ■ ■ ■ is Euler constant, and 7 ^ = as /77 is the quark anomalous dimension. 
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The h’s in the decay amplitudes are given by performing Fourier transformation on the 
transverse momenta for propagators of virtual quark and gluon in the hard part calcu¬ 
lation, they result in 


ha{x2,X3,b2,h) 


h‘'^\xi,X2,X3,bi,b2) 


{'n-i/2f Hl{MB^Jx2X3b2)St{x3) (A8) 

X [HQ{MBy/^b2)Jo{MBy/^b3)6{b2 — b^) + (62 ■*-*■ ^ 3 )} 5 
/ Ko(M^^6i), for F,>o\ 

\fRi^\MB./^ 61 ), for F, < Oy 

I^Y^o\^B\/x^bi)Jo{MBy/x^b2)9{bi - 62 ) + {bl ^ & 2 )| , 


with the variables Fi = X 2 {xi—X 3 ), F 2 = a; 2 -l-(l —a; 2 )(a;i-l-a; 3 ). And 
The threshold resummation form factor St{xi) is adopted from ref. 0 


St{x) 


21+2T(3/2 + c) 
A/?r(i -h c) 


[x(l 




Jo(^)+* Yo(^). 


(AlO) 


where the parameter c = 0.3. This function is normalized to unity. The hard scale Fs in 
the amplitudes are taken as the largest energy scale in H to diminish the higher order 
corrections: 


= max(MB^/x^, I/62,1/&3), 

(All) 

tl = max(MsA/ii, I/62,1/&3), 

(A12) 

= max{MB\/\Fi\, MBy/x2X3, 1/61,1/62), 

(A13) 

= max{MB\/^, Mb^Jx 2X3, I/61,1/62)- 

(A14) 
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